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Abstract 



We investigate the properties of different modifications to the hnear <j-model 
(including a dilaton field associated with broken scale invariance) at finite 
baryon density p and nonzero temperature T. The explicit breaking of chiral 
symmetry and the way the vector meson mass is generated are significant 
for the appearance of a phase of nearly vanishing nucleon mass besides the 
solution describing normal nuclear matter. The elimination of the abnormal 
solution prohibits the onset of a chiral phase transition but allows to lower 
the compressibility to a reasonable range. The repulsive contributions from 
the vector mesons are responsible for the wide range of stability of the normal 
phase in the {/j,, T)-plane. The abnormal solution becomes not only energet- 
ically preferable to the normal state at high temperature or density, but also 
mechanically stable due to the inclusion of dilatons. 

PACS number:12.39.F 
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I. INTRODUCTION 



Although the underlying theory of strong interactions is believed to be known, there is 
presently little hope to gain insight into the rich structure of the nonperturbative regime at 
high temperature and nonzero baryon density by solving explicitly the QCD Lagrangian. 
Presently, theoreticians try to overcome this unsatisfactory situation by pursuing mainly two 
methods: First, there is the possibility to solve QCD numerically on a discretized space-time 
lattice. Reliable results are currently available only for finite temperature and zero baryon 
density. Efforts to include dynamical fermions on the lattice are still in their infancy and de- 
mand a huge amount of computing time. The second possibility is to formulate an effective 
theory based on symmetries which hopefully reflects the basic features of QCD in a solvable 
manner. We will focus on the second approach since the consideration of symmetries and 
scaling may bring deep insight into a complex problem at low computational effort [0. 
Gell-Mann and Levy succeeded early with the second kind of ansatz, using the linear 
cr-model, in order to describe hadronic properties like pion-nucleon scattering and meson 
masses. 

For the description of nuclear matter saturation properties it is necessary to introduce vec- 
tor mesons so that the binding energy results from the cancellation of large repulsive and 
attractive contributions, in analogy to the phenomenologically successful a — cj-model [Q. 
Early attempts in that direction were done by Boguta who generated the vector meson mass 
dynamically by coupling scalar fields with vector mesons in the Lagrangian ||^. Unphysical 
bifurcations could be avoided within their approach, but one was unable to describe the 
chiral phase transition since the effective nucleon mass tended to infinity for p — > cxd. The 
solution m*^ = was confined to p = T = 0. Glendenning investigated the model at high 
temperatures and found no regime at finite density and nonzero temperatures where 
chiral symmetry is restored, because of the mechanical instability of the abnormal phase. 
Mishustin showed that one can simultaneously avoid bifurcations and describe a chiral phase 
transition at T=0, if one introduces an additional field Xy the dilaton, which simulates the 
broken scale invariance of QCD Q]. By coupling dilatons to vector mesons, one is able to 
obtain abnormal solutions a ~ without making the vector field massless. Originally, the 
dilaton field was introduced by Schechter in order to mimic the trace anomaly of QCD in 
an effective Lagrangian at tree level 

In this spirit, many authors applied chiral models to the description of nuclear matter prop- 
erties H-ini. In |12 one was even able to fit and describe finite nuclei as well as the widely 



used nonlinear version of the Walecka-model |T^. This model fails to describe the chiral 
phase transition, in contrast to 0, which exhibits a phase transition from a normal state to 
an abnormal one in the sense of Lee and Wick 



The aim of the present paper is to investigate the properties of these modified versions of the 
linear cr-model, which claim to give a satisfactory description of nuclear matter ground state 
properties, at finite temperature within the mean-field ansatz. In part || we present the 
model which incorporates broken scale and chiral symmetry. Our findings about its phase 
structure, the chiral phase transition in the (/i, T)-plane, and the temperature dependence 
of the nucleon effective mass are presented in part p| . 
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II. THEORY 



The linear cr-model introduced by Gell-Mann and Levy ||^ is extended to include an 
isoscalar vector meson uj and a scalar, isoscalar dilaton field x with positive parity. The 
scalar field a is the chiral partner of the pion and provides intermediate range attraction. 
The Lagrangian, which includes the particular ansatze of reads: 



— (^kin + (^Dirac — ^vec " Vq — VcSB 

Lkin = \d,ad^<y + ^d.Tzd^-n + ]^d,xd'x - -/..F^" (1) 



Cuirac = N{i'-i^d^ - gu.if.uj'' - ga{o- + nsTT " r)]N 
V.ec = -lu;,uj'^ml[r{-f + (1 - r){^r 



Xo 

The field strength tensor reads Ffj^^ = d^oj^ — dyOJ^. The original a- mo del is supplemented 
by nucleons which obey the Dirac equation and by vector mesons whose mass is generated 
dynamically by the a and x fields. We introduce a parameter r which allows the vector 
meson mass to be generated by a and x fields, respectively. The chirally invariant potential 
is rescaled by an appropriate power of the dilaton field x in order to be scale invariant. The 
effect of the logarithmic term ~ x^lnx is two-fold: First, it breaks scale invariance and 
leads to the proportionality 9^ ~ x^ ^is can be seen from 

.J.4r-.f-2a„.^^.^ . (2) 



which is a consequence of the definition of scale transformations |TB]. Second, the logarithm 
leads to a non-vanishing vacuum expectation value for the dilaton field resulting in sponta- 
neous chiral symmetry breaking. This connection comes from the term proportional to x^cr^- 
With the breakdown of scale invariance the resulting mass coefficient becomes negative for 
positive fco and therefore the Nambu-Goldstone mode is entered. The comparison of the 
trace anomaly of QCD with that of the effective theory allows for the identification of the 
X field with the gluon condensate: 

= (^G^^GD = (1 - 6)x' (3) 

The term ~ 5x'^ In cr contributes to the trace anomaly and is motivated by the form of the 
QCD beta function at one loop level, for details see [|12[. The last term Vcsb breaks the 



chiral symmetry explicitly and makes the pion massive. It is scaled appropriately to give a 
dimension equal to that of the quark mass term ~ mqqq of the QCD Lagrangian. 
To investigate the phase structure of nuclear matter at finite temperature we adopt the mean- 
field approximation In this approximation scheme, the fiuctuations around constant 



vacuum expectation values of the field operators are neglected: 
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a{x) = (a) + 5a^ (a) (4) 

xix) = (x) + Sx^ (x) 

The fermions are treated as quantum-mechanical one-particle operators. The derivative 
terms can be neglected and only the time-like component of the vector meson uj = (cuo) 
survives as we assume homogeneous and isotropic infinite nuclear matter. Additionally, 
parity conservation demands (tt) = 0. 

It is therefore straightforward to write down the thermodynamical potential of the grand 
canonical ensemble Q per volume y at a given temperature T and chemical potential /i: 

^ = V,ec + Vo + VcsB - V„„e " ^ / c^='A;[ln(l - Uk) + ln(l - Uk)] . (5) 

The free energy / is given by 

fiP,T;a,x,uj) = I^P+^ ■ (6) 

The vacuum energy Vvac (the potential at p = and T = 0) has been subtracted. 7 is the 
fermionic spin-isospin degeneracy factor (4 for the nuclear medium), Uk and denote the 
Fermi-Dirac distribution functions for fermions and anti-fermions, respectively: 

nk{T,fi*) = ^^^^^^^ _ ^ ^; MT, P*) = ^^^.^^^ ^ ^ ^ , (7) 



where the single particle energy is E*(k) = yk'^ + ^at^ with — QcrCr. The effective 
chemical potential reads /i* — /i — gijCU. The meson fields are determined by extremizing 

= -u;ml[r{^r + (r - 1)(^)^] + g^p ^ (8) 
d(jj (To Xo 

—-— = -uj -2 X-ko-2'^ +(4-4 + l + ln-j-251n-2)x - 2mJ,-2 = (9) 

^X Xo Xo Xo Xo "^0 Xo 

— ^ = -u; ^-(T - fco(— ) cr + Act"' - d m^U{—) + Qaps = , (10) 

9(T (7^ Xo Xo 

where the scalar density is given by 

+ ■ (11) 

The vector field cu can be solved explicitly in terms of a and x, yielding 

^ ^^(r(^)4Tr-l)(3^)^) ■ ^^^^ 

Note, that in cuo-direction the pressure is minimal, since the temporal and spatial components 
of the vector field enter with opposite sign and only the latter are dynamical variables. 
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In addition, one has to determine the baryon density at a given chemical potential via the 
equation 



P = 7 



(2^ 



The energy density and the pressure are given by 



e — Vvec + Vo + VcSB — Vmc + 
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(2vr)= 



d^kiE'^ik) - fi*) ^ esB = WsbT 



(13) 



(14) 



^ p=0 
P= ~y ' PSB 



Here, the index SB denotes the corresponding quantities in the Stefan-Boltzmann limit with 
a SB = 77r2/120. The limit T — can be taken straightforwardly, using 



lim^Tln(l - Uk) = E*{k) - fi* 



Applying the Hugenholtz-van Hove theorem [|17|, the Fermi surface is given by 

E*{kF) = ^k^p+{g^ay = fx* . 
The scalar density and the baryon density can be determined analytically, yielding 



47]-2 



kpEp — m*^ ln( 



,kF + El 



m 



N 



(15) 



(16) 



(17) 



P = 7 



''F d^k 7/c| 



(27r)3 Qtt^ 



If the dynamical vector meson mass is considered as being generated by x alone and if 6 is 
set to zero, it is possible to solve equation |10| analytically: 



X = Xo, 



\\ koa + mlf^ 



(18) 



Thus, the numerical procedure is simplified to finding the root of a nonlinear equation of 
one independent variable, namely a. This allows for a visualization of the phase structure 
at zero temperature. 

In order to describe hadrons and nuclear matter within the model, the appropriate model 
parameters must be chosen. The pion mass is fixed at the value = 138 MeV which 
determines the parameter k^ from the following relation: 



J IT 



(19) 



This equation can be obtained from equation [T0| by setting p = T = and using ctq = f-,,. 
In addition, one has to ensure that also in the vacuum ^ q-^^ = 0. This leads to the 
determination of ki: 
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h = ^{2ml + ko-p . (20) 

The Goldberger-Treiman relation can be used at the tree level to fix the coupling of the 
nucleons to the a field, = 

The vector meson mass is set to m^j = 783 MeV and A is a free parameter which determines 
the (T-mass. The remaining two parameters {g^ and xo) are fitted to the ground state nuclear 
matter binding energy Eb = e/p—rriN = —16 MeV with zero pressure at equilibrium density 
Po = 0.15 fm~^. Several parameter sets have been tested. They are listed in table 1. The 
first three rows correspond to the version of [^] with 6 = 0, which we will call hereafter 
the minimal model. There, the vector mesons are coupled only to the dilatons. Concerning 
the compressibility K, which should be around 200-400 MeV |l^, we find that a small 



quartic self-interaction of the a corresponding to small A is to be preferred in this model. If 
the logarithmic potential is included proportional to 6 |]12|, it is possible to set A = and 



therefore to lower the compressibility to reasonable values. Note, however, that the effective 
nucleon mass at po, which should be ~ 0.7mN, tends to increase with decreasing A. 



III. RESULTS 

In order to study the properties and the impact of the different modifications to the 
minimal chiral model on the observables, we focus first on the phase structure at T=0 
before discussing our findings at finite temperature. 

The infiuence of the explicit chiral symmetry breaking term on the phase structure of nuclear 
matter is checked by computing the binding energy of nuclear matter versus the a field for 
normal nuclear density p = po (Fig. 1 above) and p = Apo (Fig. 1 below) in the minimal 
version of the chiral model with 6 = 0. The first and the second column correspond to 
the model with and without explicit symmetry breaking, respectively. According to 0, the 
phase curve in Fig. la exhibits the appearance of three distinct minima: The first one is 
at m*^ ~ 0.6 — O.Tm^r (the exact value depends on the parametrization) which we denote 
as the 'normal' minimum. Besides a metastable minimum at roughly ~ 0.2m at, which 
does not play a significant role (it never becomes the energetically lowest state), there is a 
third minimum corresponding to nearly vanishing effective nucleon mass (0.02m at). This is 
the 'abnormal' minimum which becomes the energetically preferable state for large densities 
(Fig. lb). There, a phase transition takes place into a chiral phase where the nucleon 
effective mass as order parameter is nearly vanishing. Fig. Ic shows that the exclusion 
of explicit symmetry breaking effects in the Lagrangian does change the phase structure 
even at po significantly. Although the properties of the matter at the normal minimum are 
not affected, the exclusion of the explicit symmetry breaking term eliminates the abnormal 
solution entirely and therefore a chiral phase transition does not occur. 
There is another constraint for the existence of an abnormal phase: A pure cu-cr-coupling 
without a dilaton admixture (r = 1) eliminates the abnormal solution. This can be seen as 
follows: For r = 1, an additional term enters the numerator of equation IT8| yielding 



Aff3 + g^p^ - ^2p2^2^2^^3 
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so that X diverges for o" — > 0. In fact, irrespectively of which parametrization one uses, x 
becomes imaginary as soon as a ~ QAmj^. No solution is possible for smaller a values, where 
an abnormal minimum would occur. We tried to lower the compressibility in the minimal 
version of the chiral model presented in and found a lower bound of A=150 necessary to 
ensure that the abnormal state is not the energetically lowest one at normal nuclear matter 
density. If an abnormal minimum exists, at ground state density, the a^-term has to con- 
tribute strongly and the compressibility cannot be lowered to observed values. A way out is 
to permit 5 7^ |jl2[, which mimics the contribution of quark pairs to the QCD /3-function 



at one loop level: Then, it is possible to break the symmetry spontaneously even without a 
quartic self-interaction, i.e. with A = 0. The compressibility is thus lowered to reasonable 
values, without abnormal or chiral phase restoration occuring at high energy densities. 
Let us now turn to finite temperatures. Here, the analysis gets more involved: three cou- 
pled equations have to be solved simultaneously. At low temperatures, the model exhibits 
a liquid-gas phase transition as can be seen from Fig. ^ (using parameter set V). The main 
difference between the minimal and the extended model sets in at high temperatures and 
densities because of the existence of the abnormal solution in the minimal model. Fig. |^ 
shows a contour plot of the free energy at T=170 MeV and at ground state density po using 
set I. The abnormal minimum (at nearly vanishing nucleon effective mass) and a normal 
phase (at m*^ ~ 0.7m7v) are clearly visible. At normal nuclear density, a chiral phase tran- 
sition occurs at T=168 MeV. The phase transition is of first order, since the change in the 
free energy is discontinuous. 

The calculation of the phase boundary in the (/x, T)-plane yields surprising results if the 
minimal model is used (Fig. ^, Set I). Along the boundary shown in the figure the difference 
between the pressure of the abnormal and normal solutions vanishes, i.e., the transition to 
the chiral phase takes place. The transition at T=0 was already noted in [Q. 
However, the extension to finite temperatures does not lead to a closed phase boundary, 
regardless which parametrization one uses (see, e.g., triangles with A = 300, black circles 
with A = 220). The abnormal solution is stable at high temperatures or at high baryon 
densities, but not for both. This can be seen from Fig. |^, where at four particular points in 
the (/X, T)-plane of Fig. ^ the pressure as a function of the a field is drawn. The abnormal 
maximum of the pressure is flat (Fig. §a,b) or it disappears completely (Fig. ^) far away 
from the phase transition line. It becomes a well pronounced maximum with a high barrier 
to the normal state in the vicinity of the phase transition region (Fig. ^i). 
The result that one has an open phase boundary within the plotted (/i, r)-regime is unusual 
and counter- intuitive^. In contrast, in ||2^ a closed phase boundary was obtained by investi- 



gating the linear cr-model including neither repulsive contributions from cj-meson exchange 
nor dilatons. To simulate this calculation within our model we keep all parameters constant 
and change only the a;-coupling to g^i = Q (black dots) and Qui = ^ with varying gluon 
condensate (white trangles) and Quj = ^ with the gluon condensate frozen at its vacuum 
value (white circles). The presence of the dilaton field does not lead to the fan out of the 
phase transition curve. Nevertheless, it has the considerable effect to shift the transition 



^However, the increase of the critical chemical potential at small temperatures can be shown 
analytically in a low temperature expansion IH 
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points to roughly twice the values as compared to the 'non-frozen' case. Switching from 
goj = to = 6 and to g^^ = 8.2, the phase boundary spreads out to higher densities 
and temperatures. Therefore, the reason for the unusual form of the phase boundary is the 
repulsive contribution due to the cu-meson exchange. 

At that point we should emphasize that our results are obtained in the framework of the 
mean-field approximation. The inclusion of quantum fluctuations in the meson fields could 
change our findings qualitatively. This wil be investigated elsewere [^. Inclusion of reso- 



nances might lead to the closure of the boundary as was observed in ||2^ and that taking 



these additional degrees of freedom into account, the critical densities and temperatures 
decrease. Another possibility to get a closed phase boundary might be the inclusion of a 
quartic self interaction for the vector meson, (cUf^uj^y, yielding u ~ p^/^: the amount of 



repulsion at high densities is lowered. A detailed analysis will be found in . 
The extended chiral model with 6^0 does not show a chiral phase transition at all. The 
nucleon effective mass increases at high density and temperature^, as can be seen in Fig. |^. 
A similar behaviour of the effective nucleon mass can be found for the normal phase of the 
minimal model. The difference to the extended model comes from the fact that -according 
to the phase diagram of figure ^ a transition from high to low effective masses or vice versa 
can be found. 

In contrast to finite baryon density, almost no temperature dependence of the effective nu- 
cleon mass in the normal phase is found at p = until the phase transition takes place. In 
addition, the abnormal phase at p = differs qualitatively from the one at finite density. 
There, the two fields a and x vanish exactly, irrespective of the explicit symmetry breaking 
term, whereas at finite baryon density the x in the abnormal phase remains finite, as 
can be seen in Fig. |^ (there, x ~ 84 MeV). When a = 0, the scalar density vanishes and 
from equation |TB] it follows that the x fi^ld becomes zero. Because the baryonic density 
vanishes, no singularity occurs if x = 0. 

It is also interesting to compare the high temperature phase transition of the Walecka-model 
at zero density studied in with that of the minimal chiral model (Fig. One observes 



that at high temperatures the energy density and the pressure asymptotically approach the 
limit of a noninteracting fermion gas. As in we find that the energy density decreases 



with high temperatures whereas the pressure reaches its asymptotic limit from below. 
Similar results concerning the properties of the linear a-model at finite temperature were 
obtained in [^, which in our terminology would be the minimal model with a pure a — uj- 
coupling and no dilatons. However, there is an important difference, which results from the 
inclusion of the dilaton field x'- Whereas in the abnormal phase is always mechanically 
unstable (the pressure decreased with compression), leading to the result that no region in 
the (p, T)-plane existed where chiral symmetry was restored, we find here that the abnormal 
or chiral restored phase is always mechanically stable (Fig. H). The difference to Glenden- 



ning's work originates from the a; — rather than u — a-coupling. In contrast to [^], where 
it is argued that the influence of the dilaton is negligible at finite density because of its high 
mass, we find the variation of the condensate x to be essential for a mechanically stable 



^This general behavior in the chiral a — w-model is in contrast to that suggested by the Nambu- 



Jona Lasinio model [^,24| which cannot reproduce the binding energy of nuclear matter properly. 
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abnormal phase. Similar results pointing to the importance of the dilaton field in nuclear 



matter are also obtained in |2^, where the Walecka model including dilatons was studied. 
The final question to be addressed is, whether the interesting (T,p)-regions can be reached in 
relativistic heavy ion collisions. For a rough estimate, we solve the Rankine-Hugoniot-Taub 
adiabate (RHTA), which can be used as a first approximation for the description of nearly 
central collisions of fast heavy nuclei ||29|j30[| . The thermodynamic quantities calculated for 
the compression stage of the collision are shown in figure ^ The gap in the solution of the 
abnormal branch comes from the disappearance of the abnormal maximum in the corre- 
sponding region (see, i.e. Fig. ^. 

The evolution of the system in the subsequent expansion is calculated by the isentropes 
starting from a point on the Taub-adiabate (Fig. [lO|a). A minimum at S/A < 2 in the tra- 
jectory allows for the mechanical instability, which is suggested to cause multifragmentation. 
The expansion of the system from an abnormal initial state through a mixed phase into the 
normal state is shown in Fig. |TD|b. Even though we cannot reach the abnormal phase with 
the shockfront model, it might be possible, i.e. with the fireball model with p = 2^cmPo- 



IV. SUMMARY AND OUTLOOK 



The properties of the linear cr-model presented in [^,|T0|,|T5|,|T2| are studied at finite temper- 
ature T and nonzero baryon density p. At nuclear matter saturation density po, the minimal 
model of exhibits two phases (the abnormal one at nearly vanishing nucleon mass and 
the normal phase at m* ~ 0.7m7v), which allows for a phase transition at high temperatures 
or high densities. The presence of vector mesons leads to an open phase boundary, and the 
inclusion of dilatons makes the abnormal phase also mechanically stable. However, in the 
model abnormal solutions at po exist only at unphysically high values of the compressibility 
(K ~ 1400 MeV). Therefore, the abnormal phase should be eliminated by either including a 
uj — (j-coupling or by replacing the quartic self-interaction with a logarithmic term {S ^ 0). 
In this case, no chiral phase transition can be found since the nucleon effective mass as order 
parameter increases at high densities and temperatures. It remains a challenge to construct 
a reasonable chiral model for nuclear matter which allows for the study of phase transitions. 
First calculations done in an extension of the model to SU(3) are encouraging 
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FIG. 1. Binding energy versus cr/aQ with (left) and without (right) explicit symmetry breaking 
(ESB) for saturation density po (above) and 4po (below) calculated with Set I 
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FIG. 2. Liquid-gas phase transition in the chiral a — u model (calculated with Set V). 
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FIG. 3. A contour plot of the free energy in the (x, (T)-planc. The abnormal and normal 
minimum are visible. 
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FIG. 4. Phase transition points for different values of the N-w couphng constant in the (/x, 
T)-plane 
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FIG. 5. Pressure versus the scalar field a. The abnormal solution does not exist far from the 
phase transition line. The barrier between normal and abnormal phase becomes well pronounced 
near the phase transition region. 
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FIG. 6. Effective nucleon mass versus temperature (above) and density (below) 
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FIG. 7. High temperature limit of the energy density and pressure for zero density (Set I). 
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FIG. 8. Pressure as a function of density (left) and chemical potential (right), calculated with 
Set 1. 



18 




19 



600 




RHTA-Adiabate 
0=2.0, 

(7=3.0, 

(7=4.0, 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



2500 



2000 



> 1500 
<u 



< 

1000 



500 



(7=11.7 

— normal 

■ ■ ■ ■ mixed phase 

— abnormal 



0.0 



0.5 



1.0 



1.5 2.0 
p/pO 



2.5 



3.0 



3.5 



FIG. 10. Expansion of a compressed state along different isentropes starting from the normal 
phase (a) and from the abnormal phase (b), respectively. 
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TABLES 



TABLE I. Different parameter sets whieh describe nuclear matter saturation 



Set 


A 


Xo (MeV) 






X/Xo 


K (MeV) 


r 


33(5 


mT,{MeV) 


I 


300 


189.3 


8.2 


0.66 


0.71 


1464 








138 


II 


220 


188.7 


8.2 


0.67 


0.71 


1403 


0.5 





138 


III 


40 


331.7 


6.8 


0.78 


0.94 


669 


1 





138 


IV 


0.84 


392.9 


5.9 


0.84 


0.99 


387 


1 


4 





V 





372.5 


7.6 


0.80 


0.98 


356 


0.5 


4 
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